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By requiring general-coordinate and conformal invariance of the hydrodynamic equations, we show 
that the unitary Fermi gas has zero bulk viscosity, = 0, in the normal phase. In the superfluid 
phase, two of the bulk viscosities have to vanish, Ci = Ca = 0, while the third one ("3 is allowed to 
be nonzero. 
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The unitary Fermi gas is a system where fermions of 
two components interact through a zero-range two-body 
potential fine-tuned to infinite scattering length. This 
system is under active study, both experimentally [1-4] 
and theoretically. At zero temperature, the system is 
located midway between the Bose-Einstein condensation 
(BEC) and the BCS regimes [5-7]. The most remarkable 
property of this system is the absence of any intrinsic 
dimensionful parameter except the density and the tem- 
perature. This implies, on the one hand, a lack of any 
small perturbative parameter. On the other hand, the 
properties of the system are universal, i.e., independent 
of the microscopic details, which leads to great simplifi- 
cation in many problems. 

Several recent experiments have concentrated on dy- 
namic properties of the unitary Fermi gas. The trap 
breathing modes was studied in Ref. [8-11], in which it 
was observed that the damping rate reaches a minimum 
near unitary. The expansion of the gas as it is released 
from a trap was studied in Refs. [12, 13], where a hydro- 
dynamic behavior was observed. Such dynamic processes 
should depend, to one or another degree, on the kinetic 
coefficients. According to fluid mechanics [14], a nor- 
mal gas has three kinetic coefficients: the shear viscosity 
ry, the bulk viscosity C,, and the thermal conductivity k. 
This is also the situation for the Fermi gas above the 
critical temperature. Below the critical temperature, the 
Fermi gas is in the superfluid phase, and the number of 
kinetic coefficients is five [14, 15]: 77, k, and three bulk 
viscosities, normally denoted as Ci, C2, and Cs- 

The purpose of this Letter is to show that at unitarity 
(i.e., when the scattering length is infinite), certain ki- 
netic coefflcients vanish exactly. Namely, in the normal 
phase the bulk viscosity vanishes, and in the superfluid 
phase two of three bulk viscosities, ^1 and (,2, vanish. 

To emphasize the nontriviality of this result, we note 
that it is by no mean related to the well-known fact 
that the Boltzmann equation for a classical monoatomic 
gas, with two-body collision, implies zero bulk viscosity. 
When three-body collisions are included, the Boltzmann 
equation generically yields a nonzero bulk viscosity [16]. 
In our case, the bulk viscosities vanish even when the 
Boltzmann equation cannot be used, but the result is 
not expected to hold outside the unitarity regime. 



In the unitarity regime, scale invariance has been 
used to derive nontrivial relationships between thermo- 
dynamic observables [17]. But scale invariance, per se, 
does not imply the vanishing of any kinetic coefhcient. 
It only implies that all kinetic coefficients are homoge- 
neous functions of the temperature T and the chemical 
potential /z. From scale invariance one can immediately 
write down the scaling behavior of the shear viscosity ry, 
the bulk viscosity ^, and the thermal conductivity k, 

V = ^nfi [—] , C ^ hnC (—] , K = —k [ — ] , (1) 

where fj, C, and k are dimensionless function of the ratio 
of the temperature and the chemical potential. In the 
superfluid phase, instead of one function C one has three 
functions Ci, C2, Cs- 

The result advertised above is therefore stronger than 
what simple scaling arguments would imply. Hcuristi- 
cally, the vanishing of the bulk viscosities can be un- 
derstood as follow. Consider a blob of a unitary Fermi 
gas that undergoes a uniform expansion, where the ve- 
locity V at point x being v(x) = cx, where c is some 
constant. Because the unitary Fermi gas does not have 
any intrinsic scale, it should remains in thermal equilib- 
rium throughout the whole process of uniform expansion. 
This means entropy is not produced during a flow with 
V — cx. Looking at the equation for entropy production, 
one finds C = 0- Applying the same argument for a blob 
of a superfluid unitary Fermi gas undergoing uniform ex- 
pansion with the same normal and superfluid velocity 
proflle, v„ — Vs — cx, we find (2 ~ 0. From the well- 
known inequality < C2C3 concludes Ci too has to 
vanish, and only ^3 can be nonzero. 

Putting this argument on a more precise footing is, 
however, not straightforward. There is no regular solu- 
tion to the hydrodynamic equations that describe uni- 
form expansion. In this rest of this Letter we put the 
heuristic argument above on a firm mathematical ground. 
We will show how the vanishing of ^ in the normal phase 
and Ci and C2 follows from the requirement that the hy- 
drodynamic equation exhibits the conformal invariance 
of the microscopic theory. 

A general discussion of symmetries of the unitary 
Fermi gas is given in Ref. [18]; for convenience we briefly 
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review these symmetries here. Let us start by discussing 
the microscopic theory. Due to the universahty of the 
unitary Fermi gas, any short-range two-body interaction 
can be used, if it corresponds to infinite scattering length. 
In particular, we can choose to work with the following 
local Lagrangian (here and below h= 1), 

C = i^^dt^-^m' + qoi^^^a-liVaf - , (2) 
which describe a system of fermions interacting through 



the Yukawa potential V{r) 



-q^iATrry^e-''/""". Infi- 



nite scattering length is achieved by requiring mq^rQ to 
be equal to a critical number (whose numerical value is 
21.1...). Moreover and go can be time-dependent; the 
only requirement is that q^VQ is held fixed. For univer- 
sality we also need to keep 'ro(t) small compared to any 
length scale in the problem. In particular, if one perform 
the following transformation 



go(i) ^goW =7W'Zo(i), 



(3) 



then, from the point of view of the physics at length scale 
larger than ro, we have mapped our theory to itself, but 
not to another theory, provided that mq^r^ is tuned to 
the value corresponding to infinite scattering length. 

The best way to expose the symmetries of the system is 
to put it in a background gauge field and in a curved 
space with a 3D metric tensor gij{t,x) [18]. The curved 
space here is a not physical; it is simply a trick to get 
physical results. The action is now 



-I 



dtdxy/g 

+ {qoa - Ao)'^p^^|J - 



—diadja ■ 



2r2 



(4) 



where g = det \gij\ and g^^ is the inverse metric of gij. 
The action is invariant under the following infinitesimal 

transformations: (i) gauge transformations, parameter- 
ized by a function of space and time a{t, x). 



Sip = iaip, 5Ao = —dta, 



6Ai 



-dta: 



(5) 



(ii) local diffeomorphism, parameterized by three func- 
tions of space and time ^*(i, x), i = 1,2, 3, 



SAo 
5Ai 



Sa 



-edkAo-Ak^\ 
-^dkA, - Akd.e ■ 
<^dkgtj - gtkdjS!' 



'-k 



mgikS. 



(6a) 
(6b) 
(6c) 
(6d) 



and (iii) "conformal transformations," parameterized by 
a function P{t) of time only. 



SO = -PO - A[0]/30, 



(7) 



where A[0] is the dimension of the field O, defined so that 
AM = A[a] = f , A[Ao] = 1, A[A,] = 0, and A[gij] = 
— 1. One can think about (7) as time reparameterization 
t ^ t' = t + (3. By direct substitution one can check that 
the action (4) is invariant under the transformations (5) 
and (6), and under (7) combined with (3) with 7=1 — 
j$. Galilean invariance is a combination of (5) and (6) 
with a = toV • X and ^'^ = V'^t, and scale invariance is a 
combination of (6) and (7) with /3 = bt and ^'^ = ^bx'^, 
where V'^ and b are constants. 

The local diffeomorphism (6) reduces to reparameter- 
ization of space when ^'^ is time independent: the trans- 
formation laws come directly from the 3D tensor struc- 
ture of the object under consideration (i.e., scalar in the 
case of tjj, a, and Aq, vector in the case of Ai, and tensor 
in the case of g^j). Equations (6) extend this invariance 
to time- dependent diffeomorphisms. 

In flat space gtj = 5ij and zero background fields 
Aq = Ai = 0, the long-distance, long-time dynamics of 
the system is described by a set of hydrodynamic equa- 
tions [14]. This should remain true when the external 
fields is turned on and when the space is curved, if the 
following conditions are met: the fields are sufficiently 
weak and vary over length and time scales much larger 
than all microscopic scales; the curvature of space is 
small and varies slowly, also compared to all microscopic 
scales. The hydrodynamic equation should then be prop- 
erly modified to take into account the external field and 
the curved metric. In writing these equations down, we 
require that the symmetries of the microscopic theory are 
inherited by the hydrodynamic theory. This condition 
arises from the fact that the hydrodynamic equations, 
which can be used to determine the response of the sys- 
tem on external perturbations, imply concrete forms of 
the fully retarded Green's functions of the hydrodynamic 
variables. The Ward identities that come from the sym- 
metries (6) should be satisfied by these fully retarded 
Green's function, which is achieved if the hydrodynamic 
equations have the same symmetries [20]. 

Our strategy is the following. First we will modify 
the standard hydrodynamic equations for the case of a 
nonzero external gauge field and a curved metric. There 
is a unique way to do so, as we shall see. Then we 
will show that the resulting equations are inconsistent 
with the conformal invariance unless certain kinetic co- 
efficients vanish. 

Consider the normal phase first. We shall write the 
hydrodynamic equations in term of the local mass density 
p, the local velocity v^, and the local entropy per unit 
mass s. The set consists of the continuity equation, 

^5t(V5P) + "^iipv') = 0, (8) 
the equation of momentum conservation, 

■^dti^ pvi) + Vfcn«= = ^{Ei - FikV% (9) 



3 



and the entropy production equation 

■^dti^ps) + Vi{pv%s - ^O't) = ^ . (10) 

We follow closely the notations of Ref. [14]: Uik is the 
stress tensor, k is the thermal conductivity, and R is the 
dissipative function. The obvious modifications are the 
replacement of the derivatives di by the covariant deriva- 
tives Vi [19] and the appearance of the force term in the 
momentum conservation equation (9), which comes from 
the electric {Ei ~ dtAi — diAo) force and the magnetic 
Fik = diAk — dkAi) Lorentz force. The stress tensor can 
be written as 



Ilife = pViVk + PQik - (r'ik, 



(11) 



where p is the pressure and ct-^ is the viscous stress ten- 
sor. The information about the kinetic coefficients are 
contained in a'^f. and R. 

Consider the dissipationless limit first, setting a' = 
R — 0. One can check that the hydrodynamic equa- 
tions are invariant with respect to the general-coordinate 
transformations, provided that Aq, Ai, and gij transform 
as in Eqs. (6), and p, s, and transform as 



5v' = -(!'dkv'+v^dke+C- 



(12) 
(13) 



The transformation law for contains terms coming 

from the vector nature of w*, but also a term that 
can be understood if one recall v changes under Galilean 
boosts (f = VH): V ^ V V 

Now consider the dissipative terms. To keep the equa- 
tion consistent with diffeomorphism invariance, one must 
require that a[^. and R transform as a two-index tensor 
and a scalar, respectively, 

= -edko'ij - akj di^ - cJikdje , (14) 
SR = -^''dkR. (15) 

In flat space the viscous stress tensor is given by 

4 = vidiVj + djVi) + (C - |r?) SijdkvK (16) 

The naive extension to curve space, 

= V{^^V3 + ^jV,) + (C - f r?) g^,Vkv'', (17) 

is, however, not a pure two-index tensor: its variation 
under diffeomorphism contains extra terms proportional 
to ^'=: 



+ (C - h) 9^3^k^ 

+ terms in Eq. (14). 



(18) 



These terms can be canceled out by adding terms pro- 
portional to the derivatives of the metric tensor to a[f.. 



Limiting oneself to terms containing the least number of 
derivatives, the terms needed are determined uniquely. 



(19) 

Now one can check that cr-^ transforms according to 
Eq. (14). 

Similarly, the dissipative function R becomes, in 
curved space 



2-R = 2 [^i'^j + Vjt^i - ^QijVkv'^ + gij - -^Qij- 



+ ^yy + ig) +^diTd'T. (20) 

Let us now specialize on Fermi systems at infinite scat- 
tering length, and discuss its conformal invariance. The 
dissipationless hydrodynamic equation is invariant un- 
der (7), if the dimensions of different fields are 

A[p] = 2A[V] = i , A[s]=0, A[^;^] = l. (21) 

Now let us consider the dissipation terms. From dimen- 
sional analysis one find that one have to set 



A[r?] = A[C] = A[k] 



(22) 



for the hydrodynamic equation to be scale-invariant. 
However, conformal invariance is not preserved generi- 
cally. The culprit is the gij that transform as 



(23) 



which leads to cr^ and R not to conform to the pattern 
of (7), 

<54 = -Pa[j - ^alj + ICPgij , (24a) 
6R = -ak - ^pR + ^C(3 (Viv' + ^) , (24b) 

unless the bulk viscosity C vanishes. Thus the require- 
ment of conformal invariance of the hydrodynamic equa- 
tions implies C = 0. [Had we known only the scale in- 
variance, i.e., had we restricted P(t) to be of the form 
P{t) = bt, the P terms would have been absent from 
Eqs. (24) and we would have been unable to reach this 
conclusion.] 

Similarly, we can repeat the argument for the super- 
fluid case. The hydrodynamics of superfluids contains 
an additional degree of freedom, which is the condensate 
phase if, whose gauge-covariant gradient is the superfluid 
velocity. 



vt = -idi^ + Ai). 



(25) 
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It transforms in the same way as the normal velocity 
Vi = v" under general-coordinate and conformal trans- 
formations. A consequence is that the relative veloc- 
ity between the siipcrfluid and the normal component 
ui' = — transforms as a pure vector under diffeo- 
morphism, 

5w' = -i^dkw' + w^dkC- (26) 

The term in the variation cancels between 5vs and 5v. 

The diffeomorphism-invariant dissipative function in 
curved space is 

+ 2Cl {WiV' + ^) VjiPsW^) + C2 {^^V' + ^) 

+ (:3[yi{Psw')? +'^diTd'T. (27) 
Under conformal transformations, R transforms as 

5R = ^Ci/3Vi(p«;0 + ^C2/3 (viv' + ^ j . 

(28) 

The requirement of conformal invariance of superfluid hy- 
drodynamics implies that the (3 terms must have vanish- 
ing coefficients, i.e., ^1 = C2 = 0. 

In conclusion, we find that in the unitary limit the bulk 
viscosity vanishes in the normal phase. In the superfluid 
phase two of the three bulk viscosities vanishes. This 
vanishing of the bulk viscosities are directly related to 
the conformal invariance of the unitary Fermi gas. 

It should be possible to check the result derived in 
this Letter by using the Boltzmann equation in the two 
regimes where it applies: the high-tcmpcraturc regime 
T ^ jji and the low temperature regime T <C /U. In the in- 
termediate regime T fj, the Boltzmann equation is not 
reliable, since wc do not have weakly coupled quasipar- 
ticles. The result derived in this paper, however, should 
be valid for all regimes of 

With respect to the recent experimental findings [8-11] 
it is tempting to speculate that the dip in the damping 
rate of the radial breathing modes near the Feshbach res- 
onance could partially be due to the vanishing of the bulk 



viscosities at unitarity. Another source for the reduction 
of damping is probably the decrease of the shear viscos- 
ity and the thermal conductivity at strong coupling. A 
study of the breathing modes using two-fluid dissipative 
hydrodynamics may enable one to extract kinetic coeffi- 
cients from experimental data, and ultimately verify the 
results derived in this Letter. 
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